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In our earlier work [M. Hafez, et al, Phys. Lett. A 373 (2009) 4479] we employed the flow equation method 
to obtain a classical effective model from a quantum mechanical parent Hamiltonian called, the ionic Hubbard 
model (IHM). The classical ionic Hubbard model (CIHM) obtained in this way contains solely Fermionic oc- 
cupation numbers of two species corresponding to particles with j and J, spin, respectively. In this paper, we 
employ the transfer matrix method to analytically solve the CIHM at finite temperature in one dimension. In 
the limit of zero temperature, we find two insulating phases at large and small Coulomb interaction strength, U, 
mediated with a gap-less metallic phase, resulting in two continuous metal-insulator transitions. Our results are 
further supported with Monte Carlo simulations. 

PACS numbers: 71.30.+h, 68.35.Rh 



INTRODUCTION 

To understand the magnetism and magnetic phenomena, 
one of the basic interactions is the exchange mechanism which 
is deeply rooted in Coulomb interactions and quantum me- 
chanical indistinguishability. Therefore a fair understanding 
of the magnetic behavior of materials is not possible without 
investigating appropriate quantum spin models. Introducing 
uni-axial anisotropy to the Heisenberg model amounts to sup- 
pression of transverse quantum fluctuations (S + S~ +S~S + ), 
leading to the so called Ising model The resulting Ising 
Hamiltonian turns out to contain the basic magnetic phases 
of the original Heisenberg model, namely ferromagnetism, 
and antiferromagnetism 12J]. Although many of the interesting 
possible aspects such as spin liquid phases, spin-wave excita- 
tions, etc. yfl can not be captured by the Ising model. Classi- 
cal Ising model has the merit of being much simpler to solve, 
and admits analytical (HI 01 and graphical solutions J2l in 
various geometries JHt] in one and two dimensions which is 
lacking in the original quantum Heisenberg model. 

Dielectric properties are among the most important proper- 
ties characterizing materials. The question can be asked here, 
is there any Ising-like model that can provide basic informa- 
tions about their phase diagram, and at the same time being 
simple enough to allow for analytical solutions? We have 
taken the example of the (IHM) latZl]. This model was intro- 
duced to study the neutral-to-ionic transition in organic com- 
pounds, as well as, understanding the role of strong electronic 
correlations in the dielectric properties of metal oxides flUtl- 
This model is as follows: 
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where Cj CT (cL) is the usual annihilation (creation) operator 
at site i with spin a. U is the on-site Coulomb interaction 
parameter, and A denotes a one-body staggered ionic poten- 
tial. The kinetic energy scale is given by the real hopping 
amplitude t which prefers to gain kinetic energy by spreading 
the wave-function over the whole system, leading to quan- 
tum fluctuation of the charge density. The zero temperature 
phase diagram of this model contains Mott and band insulat- 
ing states when the energy scales corresponding to U or A 
dominate, respectively QSQ . When these two scales become of 
the same order of magnitude, the nature of intermediate phase 
still remains controversial. Some authors argue that the inter- 
mediate state is a spontaneously broken symmetry phase 0], 
while some others argue that the phase in between is metal- 
lic Q 0, OH 1J.]. In our previous investigation |@1 we em- 
ployed the method of flow equations for the quantum Hamil- 
tonian (dJ to obtain an effective Hamiltonian in which the hop- 
ping t term is renormalized to zero, producing a new nearest 
neighbor Coulomb attraction. The resulting classical Hamil- 
tonian is of the following lattice gas form: 
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where the renormalized parameters A, U, V have closed form 
expressions in terms of the bare parameters U, A fH]. We take 
the kinetic energy scale t as the unit of energy. 

By calculating the spin and charge gaps, we showed that 
this simple classical model is capable of capturing the physics 
of a metallic state sandwiched between two distinct insulat- 
ing phases at zero temperature as one increases U for a fixed 
value of A HI . Here n ia contains two species (or colors) cor- 
responding to <T = j, 1, respectively. Since n a for each "color" 
a is either or 1, it is an Ising-like variable. Therefore our 
effective CIHM Hamiltonian can be thought of, as two inter- 
penetrating Ising models on a lattice. Being a natural exten- 
sion of two species lattice gas model, allows for analytical 
solution in one spatial dimensions (ID). Here we can employ 
the transfer matrix method, constructed in terms of 3 x 3 ma- 
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trices to calculate the thermodynamic properties of this model 
at finite temperatures. At high temperatures where the ther- 
mal fluctuations wash out the quantum effects, we expect the 
results of our calculations to be a good description of the class 
of materials modeled in terms of 1HM. Also a novel graphi- 
cal solution in 2D similar to Feynman's construction can be 
worked out 112 ] . 

The paper is organized as follows: In the next section we 
calculate the grand canonical potential for this model in one 
dimension and discuss the particle density and ionicity in var- 
ious fillings. It is followed by a section focused on half-filling 
situation, and calculate the specific heat, compressibility and 
ionicity to assess the nature of the phases in the parameter 
space of U, T. Throughout the paper we have fixed the value 
of A = 20 JH]. The final section is devoted to summary and 
discussion. 



GRAND POTENTIAL 

In this section we calculate the grand canonical potential 
that can be derived from the grand partition function (GPF) 
which is defined as follows: 
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where T is Temperature, L is lattice size, and \i is chemi- 
cal potential. The {ma} denotes all possible configurations 
of the occupation numbers, which must be summed over. 
N({n,icr}) = J2ia n i° * s ^ number of particle and E ({ma}) 
is the energy that is defined by the model Hamiltonian Eq. 0. 

The summation in Eq. (0 can be calculated analytically as 
follows. The values of m„ are only zero and one, which is 
the Fermionic memory left in the commuting variables rii a . 
Hence nf a = m a , and the second term in Eq. (01 becomes: 



*T n a = \ + ™ii)[( n *T +riii)-l]. 
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The GPF becomes, 
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where the coefficient 17^=1 (•"■ + takes into account the 

two-fold degeneracy for m = 1 in Eq. (|5), which corresponds 
to either = 1, na = 0, or n,t = 0, mi = 1. Eq. (|6]i can 
be written in a more symmetric form, 
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where L is assumed to be even and the periodic boundary 
conditions, til+i = n\ is implied. Defining the matrix ele- 
ments of the transfer matrix M ni) „ 2 by, 
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Eq. <J2J can be written as: 
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where A/' is the transpose of M, M„ 2iJll = (M') nij „ 2 . Then 
the partition becomes: 



Z = Aj 2 + A 2 2 + A 
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where Ai, A2, and A3 are the eigenvalues of M t M. In the 
thermodynamic limit where L is very large, the grand poten- 
tial per site becomes: 



As can be seen at the Hamiltonian level, only the spin- 
summed occupations appear and therefore this classical model 
does not contain spin-polarized (magnetic) solutions. So we 
change the summation variable from n,f and n,j to m = 
n. L \ + mi that has three possible values 0,1,2. Therefore we 
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Here X m ax is the maximum eigenvalue. Hence to obtain 
the grand potential one needs to calculate the eigenvalues of 
M'M. Later on, when we discuss the ionicity, we will in 
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FIG. 1: Particle density n, versus chemical potential /i for U = 
10, A = 20 and different temperatures. Three plateaus at n — 
0.5, 1, 1.5 correspond to quarter-,half- and three quarter-fillings. For 
lower temperature, the variations of the density as a function of \i is 
very slow over the plateaus, which is clear signature of a gap in the 
spectrum at these commensurate fillings. By increasing the temper- 
ature, plateaus get rounder, which indicate the thermal energy starts 
to overcome the gaps. The inset magnifies the "isosbestic" behavior 
at n = 1 plateau. 



addition need the eigenvectors of M*M. The eigenvalues of 
M*M are solutions to the third order equation, 



\ A + a 2 X + aiX + a = 0, 
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with 0,2, ai, and an given in the appendix. Through the pa- 
per we will report the numerical plots for A = 20 in units 
in which t = 1. For this value of A, there would be no 
level-crossing among the eigenvalues A of the matrix M'M 
when one varies \i, T, and U, as expected from Perron's theo- 
rem 1 1311 . 

From the partition function, one can in principle calcu- 
late various averages of the form (n iry nj a i ...). The sim- 
plest of these averages are the average particle density n = 
( n 2i+i + n 2i)/% (symmetric combination), and the ionicity 
/ = (n,2i+i — «2i)/2 (antisymmetric combination) which 
contain important information about the nature of the thermo- 
dynamic phases of the model. In the following we calculate 
these averages as a function of \i. 



Particle Density 

Once the grand potential <fi is known, one can derive many 
other thermodynamic quantities. The particle density, n, can 
be calculated as, 



(13) 
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In Fig.Q] the particle density versus [i has been plotted for 
different values of T. In this figure the value of U is fixed 



FIG. 2: Particle density n, versus chemical potential [i for constant 
T = 0.4, A = 20 and different values of U. As can be seen in gen- 
eral again there are three plateaus corresponding ton = 0.5, 1, 1.5. 
However, for intermediate values of U, the plateau at half-filling has 
disappeared. 



to be U = 10. This value of U corresponds to a band insu- 
lating regime at zero temperature jjj. As can be seen in the 
figure, apart from trivial cases of the empty (n = 0) and the 
fully-filled {n = 2) lattice, there are three plateaus. The first 
one corresponds to the half-filling n = 1 (i.e. one particle per 
lattice site), and two others correspond to commensurate fill- 
ings, n — 0.5 (quarter-filling), n = 1.5 (three-quarter-filling) 
which is quite similar to the plateaus of the parent ionic Hub- 
bard model floll . Indeed the n = 0.5 and n = 1.5 are related 
by a particle-hole transformation. Therefore in what follows, 
we focus on n = 0.5, 1. With increasing the temperature, the 
plateaus get more and more rounded. The value of chemi- 
cal potential corresponding to the half-filling is independent 
of temperature, as can be seen in Fig. Q] The inset plot indi- 
cates this point more transparently. This "isosbestic" behavior 
is observed for all values of values of 10 < U < 40. 

In Fig. |2] we have plotted the particle density as a func- 
tion of chemical potential /i at a constant (low) temperature 
T = 0.4, for different values of U. As can be seen in the 
figure, increasing U, causes the n = 0.5, 1.5 plateaus which 
are particle-hole counterpart of each other, get wider. How- 
ever the plateau at n = 1.0 gets narrower and finally vanishes 
around U = 20, resulting in a gap-less phase at half-filling. 
Upon further increasing U, the half-filling plateau is recov- 
ered, and gets wider, indicating the emergence of a growing 
new gap in the system, which is reminiscent of the Mott insu- 
lating behavior. To identify the nature of gap at n = 0.5, 1.5 
one needs to calculate the ionicity, which provides informa- 
tion about how the unit cell is being filled. This will be done 
in the sequel. 

For later reference, we report the value of chemical poten- 
tial corresponding to half-filling, which by examining the nu- 
merical plots turns out to be: 



H{T, U,n=l) = 0.495J7. 
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FIG. 3: The ionicity per site as a function of chemical potential 
for various temperatures at U — 10. Three plateaus around \i « 
-5, 5, 15 correspond to 1/4,1/2,3/4 filling. 
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FIG. 4: Ionicity per site as a function of fj, and constant temperature 
T — 0.4, for various values of U. Three plateaus correspond to 
1/4,1/2,3/4 filling, respectively. Note that the center of the half-filling 
plateau is at fi — U/2. 



The temperature range at which the above relation is valid, is 
roughly 0.1 < T. For temperatures outside this range, there 
will be deviations from the above rule of thumb relation. 



Ionicity 

If A, B denote the sublattice of odd and even sites, respec- 
tively, then Na and Nb are defined as: 
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the ionicity per site becomes I = (Na — Nb)/L, where rii = 
riij + rij|, and L is the total number of lattice sites. This 
quantity can also be calculated analytically with the aid of 
transfer matrices. In the grand canonical ensemble (rij) is 
given by: 

= I y n . e -K B (^»-^({^})). (i7) 

Depending on whether j is odd or even, similar to calculations 
for grand partition function, we obtain: 

(n 3 ■) = -^Tr (V(M*M)*) , foroddj (18) 

= -Tr(iV(MM')^, forevenj (19) 



where TV is the following 3x3 matrix: 

-Nni,ri2 — ^l^ni,n2- 
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In Fig. [3] we have plotted the ionicity per site as a function 
of fi for various temperatures at a fixed value of U = 10. 



For this value of U, at zero temperature one would expect 
a band insulator at half-filling, where A-sublattice (with on- 
site energy — A/2) are doubly occupied. Let us focus around 
\x = — 5 which corresponds to quarter filling (c.f. Fig. Q]). 
The value of I = 0.5 for /i w —5 in Fig. [3] shows that the 
only electron of each unit cell belongs to A-sublattice, which 
represents a charge density wave insulator. 

As we increase /i from « — 5, the ionicity increases, which 
indicates that sublattice A continues to be filled. When /j, sw 5, 
one approaches the half-filling (Fig. Q]), where at lower tem- 
peratures the A-sublattice is totally occupied, hence I w 1. 
As can be seen in Fig. [3] this saturation value is decreased as 
the temperature is increased. This reduction in the ionicity, 
indicates that the particles are bing thermally excited across 
the band gap. Further increasing of the chemical potential, 
one reaches the plateau around // w 15 of Fig. [3] (c.f. Fig.Q]). 
The decrease in the ionicity indicates that the added particles 
essentially start to occupy the B-sublattice. The symmetry of 
Fig.[3]around n = 1, is due to the apparent particle-hole sym- 
metry of the Hamiltonian. 

In Fig.[4]we plot the ionicity per site at the constant temper- 
ature T = 0.4 and different values of U indicated in the leg- 
end. Again there are three plateaus corresponding to 1/4,1/2 
and 3/4-filling, respectively. This figure indicates that the cen- 
ter of the half-filling plateau is at = U /2. This can be 
understood as a Hartree like energy for the Hubbard model. 
Now concentrating around half-filling in this figure, for small 
values of U, the ionicity reaches 1, which indicates a com- 
plete electric polarization of the unit cell, and hence we have 
a band insulator. For intermediate values of U, the maximum 
value of ionicity does not reach 1, which shows some of the 
added particles start to occupy the B-sublattice. Also note that 
for the intermediate values of U the half-filling plateau starts 
to disappear, i.e. the emergence of a metallic behavior. For 
large values of U, the half-filling plateau is restored, but at 
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FIG. 5: tonicity and specific heat per site at half-filling from trans- 
fer matrix and Monte Carlo (MC) methods. To stay in the range of 
validity of Eq. d 14t . transfer matrix results are plotted for T > 0.1. 
Monte Carlo results are obtained for L — 500 sites. As stated earlier, 
we are working at constant A = 20. 



FIG. 6: Specific heat per site versus T at A = 20 for different values 
of U. For small values of U there is a single peak, which becomes 
a peak-dip-hump in intermediate U, and finally the peak merges into 
the hump at large U. For U = 20, both analytic and Monte Carlo 
results are reported . MC results are for L — 500 sites. 



zero ionicity. Therefore we have an insulator with unpolar- 
ized unit cell. Such an insulating state can be thought of as 
classical counterpart of the Mott insulating state. 



HALF-FILLING 

As we saw in Fig. [2] for n = 0.5, 1.5 the U, A energy scales 
cooperate with each other, to give rise to a charge density wave 
insulating ground state. Rather at n = 1.0, these two energy 
scales compete against each other to destroy the insulating be- 
havior for intermediate U (~ A), giving rise to a richer phase 
diagram. Therefore in this section we focus at half-filling and 
calculate the specific heat, compressibility, and ionicity. Be- 
fore doing so, we compare some physical quantities evaluated 
by a fully numerical Monte Carlo simulation, with our exact 
transfer matrix results. In Fig. [5] ionicity and specific heat per 
site are plotted at half-filling and show a good agreement be- 
tween analytical and numerical results. This ensures that both 
Monte Carlo and analytic results are reliable. 

Now let us proceed with the calculation of various thermo- 
dynamic quantities. For the fixed value of A — 20, we have 
two ways to plot thermodynamic quantities. First way is to 
plot them as a function of temperature T, at some selected 
values of U. These results indicate that in the present one 
dimensional model, there will be no finite temperature phase 
transition. This is obviously due to the analytic behavior of 
the partition function as a function of T. The second way is 
to plot them as a function of U, for some selected tempera- 
tures. This second way of presenting the data, reveals that as 



one lowers the temperature, there will be sharper features as 
a function of U, indicating the zero temperature phase transi- 
tion. Our calculations are based on Eqs. dl 11141 ). For very low 
temperatures, where the validity of Eq. ( fl4l ), might be ques- 
tioned, we employ Monte Carlo simulation data. 



Specific heat 

Specific heat per site can be calculated as: 



where s is entropy per site at half-filling that can be derived 
from the grand potential. Figure [6] shows the specific heat 
versus T for various values of U. As can be seen in Fig. [6] 
for values of U < 14, there is a single peak in the cj,. For 
14 < U < 38 a peak-dip-hump structure can be observed. 
For U > 20 the hump is quite clear, while for U < 20, it can 
be interpreted as a precursor to the hump. For 38 < U, the 
peak merges into the hump structure. In terms of the parent 
quantum Hamiltonian, such hump structure indicates incoher- 
ent excitations. For U = 20, Eq. (fl4l is not reliable at very 
low temperatures. Therefore we report Monte Carlo simu- 
lation results which indicates a very sharp peak at T « 0. 
As U moves towards w 20 from both sides (lower and upper 
panels), the peak gets sharper and moves towards lower tem- 
peratures. This indicates that in T — > limit one expects a 
transition between gaped and gap-less states. This behavior 
can be simply understood in terms of a two-state model with 
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FIG. 7: Estimate of gap using two-state model as a function of U at 
A = 20. 



level spacing 6, whose specific heat is given by 



cl = k B ( x 



where x 



{l + c x ) 2 J ' 
8/{k B T). Behavior of Eq 



(22) 



for x ^> 

,2 



1 is 

like ~ x 2 e~ x , while for x <C 1 it vanishes as ~ x z . For 
the intermediate region a Schottky peak around x poa k ~ 1 
(6 ~ fcflT) arises in the specific heat. Fitting the specific 
data to Eq. d22l . in Fig. [7] we have plotted the estimated gap 
versus {/, which indicates two gaped phases. According to the 
above two-state formula, this peak corresponds to S/ksT ~ 
1, from which a gap of <5 ~ 10 -2 can be estimated. If we 
extrapolate the estimated gap for 18 < U < 22, zero gap 
region is expected to occur for U Cl = 19.30 < U < U C2 = 
19.75 which is compatible with our previous work |8[]. 

Now let us look at the specific heat data from a different 
angle. In Fig.|8]we plot the specific heat versus U for differ- 
ent values of temperatures. As can be seen, there are two peak 
structures at all temperatures, which get sharper and by lower- 
ing the temperature, they tend to accumulate around U ~ 20. 




FIG. 8: Specific heat per site versus U for different values of T. 
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FIG. 9: Compressibility versus T for different values of U. 



Extrapolating the trend of this double-peak structure to the 
limit of T — » 0, suggests two phase transitions at U c \ and 
U C 2 Si, compatible with the behavior of the vanishing gap 
region in Fig. [7] The characteristic quadratic behavior around 
U « 20 seen in Fig.[8j which according to the two-state model 
is expected to be like cl oc (S/T) 2 , indicates two continuous 
metal-insulator transitions, with 6 ~ \U — U C A, i = 1, 2. 



Compressibility 

Another quantity that can be treated in our consideration 
CIHM is the compressibility which is given by: 



1 dn 
n 2 dfi 



(23) 



where n = N / L is the density of particles per site. 

In Fig. [9] we plot the compressibility as a function of T for 
different values of U. Zero compressibility is a characteristic 
of gaped states. As can be seen in this figure for small value 
of U the range of temperatures at which the compressibility 
is close to zero is substantial, which means that the gap is so 
large that up to such temperature the insulating behavior is still 
manifest. By increasing U, this temperature range shrinks and 
becomes smaller and smaller, until around U = 20, it extrap- 
olates to zero. Increasing U beyond 20, again recovers a finite 
temperature range in which the compressibility is zero. This 
behavior confirms that a gap-less state is sandwiched between 
two gaped states. 

To see the above statement more clearly, in Fig.[l0]we plot 
compressibility as a function of U for selected temperatures. 
As can be seen in this figure, there is a region with non-zero 
compressibility, which characterizes a gap-less phase. Out- 
side this region, the compressibility decays to zero. By de- 
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FIG. 10: Compressibility versus U for different values of tempera- 
ture. 



creasing T, the width of the compressibility peak becomes 
smaller and while the height of the peak diverges as T — > 0; 
a typical characteristic of a continuous phase transition. This 
confirms the existence of a metallic phase at zero tempera- 
ture J3|. The effect of thermal fluctuations is to smear the 
edges of metallic region. This is quite intuitive, as for values 
of U near the zero temperature boundary of metallic phase 
with neighboring insulating phases, the gaps are small, and 
hence the thermal energy can overcome the gap. 

Ionicity 

In Fig. QT| we plot the ionicity per site for A = 20 and 
various values of U at half-filling. As can be seen in the figure, 
by lowering the temperature, for U < 18, the ionicity tends 
to I = 1; for U > 20 it reaches a zero temperature value of 
1 = 0; while for U « 20, it reaches a value between these two 
limits; characterizing a phase with charge fluctuations. 
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FIG. 12: Ionicity versus U for some values of temperature at half- 
filling and A = 20. 



This behavior can be understood as follows: At small U 
regime the unit cell is fully polarized at low temperatures, 
with both | and J, particles occupying the A-sublattice. As the 
temperature is increased, some of the particles get excited to 
B-sublattice by absorbing the thermal energy, ksT . Similarly 
for large values of U, at lower temperatures the unit cell is not 
polarized, due to the Coulomb term U. By increasing the tem- 
perature, thermal excitations with doubly occupied sites will 
be created, thereby increasing the ionicity. For intermediate 
U « 20, the weights of polarized and unpolarized configu- 
rations in the unit cell become comparable; hence giving the 
ionicity < I < 1. 

In Fig. [12] we plot the ionicity at half-filling as a function 
of U for A = 20 and various values of T. As can be seen 
at all temperatures, the ionicity smoothly varies between 1 for 
small values U, and for large values of U. The width of the 
transition region decreases by lowering the temperature, and 
is expected to approach the figure 2 of Ref.0. 
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FIG. 1 1 : Ionicity versus T for A = 20 and different values of U at 
half-filling. 



CONCLUSION 

In this work we studied a classical model consisting of two 
Ising like variables on a one dimensional chain. Despite the 
simplicity which essentially results form the lack of Fermionic 
minus sign problem, our model captures some of the interest- 
ing properties of the ionic Hubbard model. Various thermody- 
namic properties, such as, specific heat, ionicity, particle den- 
sity and compressibility when viewed as a function of U in a 
given temperature, indicate the presence of two gaped states at 
small and large values of U, with a gap-less state sandwiched 
between them (around U /A « 1). When the same quanti- 
ties viewed as a function of temperature, there is no sign of 
phase transition down to zero temperature. The three phase 
scenario of the zero temperature with clear zero temperature 
phase transition boundaries at U Cl (A) and U C2 (A) extrapo- 
lates to higher temperatures. However, the boundaries get 
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smeared due to thermal fluctuations, giving rise to a cross- 
over behavior. 

Mapping of D dimensional quantum models to D + 1 di- 
mensional classical Hamiltonians is a well known paradigm 
in statistical physics. Our flow equation approach [8] suggests 
an alternative approach to construct "D" dimensional classi- 
cal models which might be useful in capturing basic aspects 
of the original quantum Hamiltonian. 
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APPENDIX 
Coefficients of Eq. CD 

The coefficients ai, and ao of Eq. (fT2l i are as follows: 

a 2 = - {c~ p ( 2 A+2 ^ +U ) + 2 e~ 3/2 ( A+2 ") + cT p ( 4 

+ 2 e -V2 (3 A+2 ii+2 U+8 V) , ^ g -/3 (A+2 p+2 V) 

+ 2 e" 1/2 13 (^ +6 p ) + e~ p ( 2 ^ +8 + e' 13 (^ +2 M ) 
+ 2e" 1/2/3 ( 2 ^ +8 ^ +2Al+i )}e ,3 ( 4A1+ ^), (24) 

ai = _{_4 e -/3(2A+2 AI +c7+2 V) + 4c -l/2/3(4C/+12y+2 A1 +A) 
_ 4e -/3(2C/+10 V+A) _ 2e -l/2/3(4C7+16V+3A+2 Al ) 
_ g -/3 (2 U+8 V+2 /i+A) _|_ ^ g-1/2 (2 [7+4 V+A+6 
_ g -/3 (2 U+A+2 p.) , ^ e -l/2 /3 (4 [7+12 V+3 A+2 
1 4e -l/2^(2C/+3A+6/i+4y) _ 2 c -1/2/3(A+6m+2(7) 
+ g & -0 (2 A+2 p+U+3 \?) + g c -/3 (4 ju+A+V) _ 4 e -/3 (4 M +A) 

- 4 e"' 9 ( 2 ^ +8 ^ + ^) - 2 c~ 1/2 13 ( 4 ^ +8 ^ +2 M ) 

, g e -/3 (2 [7+9 i>+A) _ 2 e -l/2 /3 (4 !7+8 i>+A+2 M ) 

+ 2 e _|3 ( 2 ^+^+ 2 ^) - 4 c~ /3 ( 2 ^ +2 f J -+^+ 4 v) 

_ 4e -/3([7+2 V+2fj.) _ 2 c -l/2 0(2U+3A+6 fJ .+8V) 

_ 2c -l/2^(4[7+16V > +2^+A) _ 2 c -1/2/3(2[7+6m+8V > +A) 

_ 2e -l/2f)(2i/+3A+6 f .) _ 4c -/3(4 M +A+2F) 

+ 8 c^ 13 ( 0+3 ^+ 2 M ) - 4 e"' 3 (^ +4 ^+ 2 }c^ ( A+6 M ) , (25) 



a = 4e 
+ 10 c 



-2/3([/+V-3/i 



-4/3 V 



4e 



){-4. 

-/3v 



-2/3 V 



4e 



-5/3 V 



4e 



4e" 



-7/3 V 



4e 



-6/3 V 



-8/3 V 



!}• 



-3/3 V 



(26) 
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